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ABSTRACT 

We study the influence of the propagation effects on the mean profiles of radio 
pulsars using the method of the wave propagation in the inhomogeneous media de- 
scribing by Kravtsov & Orlov (1990). This approach allows us firstly to include into 
consideration the transition from geometrical optics to vacuum propagation, the cy- 
clotron absorption, and the wave refraction simultaneously. In addition, non-dipole 
magnetic field configuration, drift motion of plasma particles, and their realistic en- 
ergy distribution are taken into account. It is confirmed that for ordinary pulsars 
(period P ~ 1 s, surface magnetic field Bq ~ 10^^ G) and typical plasma generation 
near magnetic poles (the multiplicity parameter A = ne/ncj ~ 10^) the polarization 
is formed inside the light cylinder at the distance Tq^c ^ lOOOi? from the neutron 
star, the circular polarization being 5-20% which is just observed. The one-to-one 
correspondence between the signs of circular polarization and position angle {p. a.) 
derivative along the profile for both ordinary and extraordinary waves is predicted. 
Using numerical integration we now can model the mean profiles of radio pulsars. It 
is shown that the standard /S-shape form of the p. a. swing can be realized for small 
enough multiplicity A and large enough bulk Lorentz factor 7 only. It is also shown 
that the value of p. a. maximum derivative, that is often used for determination the 
angle between magnetic dipole and rotation axis, depends on the plasma parameters 
and could differ from the rotation vector model (RVM) prediction. 
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1 INTRODUCTION 

More than forty years after discovery, our understand- 
ing of pulsar phenomenon leaves an ambiguous impression. 
On the one hand, the key properties were understood al- 
most immediately (see, e.g., the monographs by Manchester 
& Taylor 1977; Lyne & Graham-Smith 1998): the stable pat- 
tern of radio emission is related to the neutron star rotation, 
the energy source is the kinetic energy of its rotation, and 
the release mechanism has the electromagnetic nature. The 
mean profiles of radio pulsars are well described by the hol- 
low cone model (see below). Within four last decades enor- 
mous amount of observational data concerning polarization 
and other morphological properties of mean profiles was col- 
lected (Rankin 1983, 1990; Johnston et al. 2007; Weltevrede 
& Johnston 2008; Hankins & Rankin 2010; Keith et al. 2010; 
Yan et al. 2011). However, there is no agreement on the 
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mechanism of coherent radio emission (Beskin 1999; Usov 
2006; Lyubarsky 2008). 

Any self-consistent theory of pulsar radio emission must 
include at least three main elements. First, it should describe 
a plasma instability that produces coherent radio emission. 
Second, the saturation of this instability that determines the 
intensity of the outgoing radio emission. Third, once the ra- 
diation is produced, its polarization properties are modified 
by the interaction with magnetospheric plasma and fields. 
These propagation effects should be accounted for to make 
a quantitative comparison of the theoretical predictions for 
the generation of radio emission with observational data. 

There are several proposed mechanisms for the initial 
instability: an unstable flow of relativistic electron-positron 
plasma flowing along curved magnetic field lines (Goldre- 
ich & Keeley 1971; Blandford 1975; Asseo et al. 1980; Be- 
skin, Gurevich & Istomin 1993, hereafter BGI); an instabil- 
ity caused by boundedness of the region of open field lines 
(Luo, Melrose & Machabeli 1994; Asseo 1995); an instabil- 
ity connected with kinetic effects, that can be caused by 
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non-equilibrium of the particles energy distribution function 
(mainly, anomalous Doppler effect in the region of cyclotron 
resonance); two-stream instability (Kazbegi, Machabeli & 
Melikidze 1991); the instability connected with the nonsta- 
tionarity of plasma particle production in the region of its 
generation (Lyubarskii 1996). 

The saturation mechanism, whose investigation requires 
involving the effects of a nonlinear wave interaction, is the 
most complex from the theoretical point of view. There- 
fore, it is not surprising that only a few researchers have 
managed to consider this question consistently (see, e.g., Is- 
tomin 1988). Finally, the processes of the wave propagation 
in pulsar magnetosphere have not yet been investigated with 
sufficient detail either, although the part of theory that in- 
cludes the propagation processes can be constructed using 
the standard linear methods of plasma physics. 

Recall that the present interpretation of the mean pro- 
files of radio pulsars is based on the ground of so-called 
hollow cone model (see, e.g., Manchester & Taylor 1977). 
Within this approach the directivity pattern is assumed to 
repeat the profile of the number density of secondary parti- 
cles outflowing along the open field lines. As the secondary 
plasma cannot be generated near the magnetic pole (where 
the curvature photons radiating by primary beam propagate 
almost along magnetic field lines), the particle number den- 
sity is to have the 'hole' in its space distribution (Sturrock 
1971; Ruderman & Sutherland 1975). 

There are four assumptions in the hollow cone model: 
first, the emission is generated in the inner magnetospheric 
regions (where the magnetic field may be considered as a 
dipole); second, the emission propagates along the straight 
line; third, the cyclotron absorption may be neglected; and 
fourth, the polarization is determined at the emission point. 
Such basic characteristics of the received radio emission al- 
low to determine the change of the position angle (p.a.) of 
the linear polarization along the mean profile (Radhakrish- 
nan k Cocke 1969) 



p.a. = arctan 



sin acosCcos(/) — sin ("cosa 



(1) 



Here a is the inclination angle of the magnetic dipole to the 
rotation axis, C is the angle between the rotation axis and 
the observer's direction, and is the pulse phase. 

As a result, the radiation beam width Wr itself and 
its statistical dependence on the period P can be quali- 
tatively explained under these assumptions (Rankin 1983, 
1990). As the pulsar radio emission is highly polarized (Lyne 
& Graham-Smith 1998), one could check the validity of the 
relation (IT]) as well. As is well-known, in many cases the 
observed p.a. swing is in good agreement with this theoret- 
ical prediction. Besides, the polarization observations show 
that the pulsar radio emission consists of two orthogonal 
modes, i.e., of two components which p.a. differ by 90° 
(Taylor & Stinebring 1986). It is logically to connect two 
such components with two normal modes, ordinary and ex- 
traordinary ones, propagating in a magnet oactive plasma 
(Ginzburg 1961). It is not surprising that the hollow cone 
model in its simplest realization is currently widely used 
for quantitative determination of the parameters of neutron 
stars. 

At the same time, it is well known that, in general, 
three main assumptions are incorrect. First of all, after the 



paper by Barnard & Arons (1986), it became clear that the 
ordinary wave (i.e., the wave which electric field belongs to 
the plane containing the wave vector k and the external 
magnetic field B) does not propagate in a straight line, but 
deflects away from the magnetic axis. Subsequently, this ef- 
fect was studied in detail by Petrova & Lyubarskii (1998, 
2000), it was an important element in the BGI theory. The 
correction to relation (fl]) connected with the aberration was 
determined by Blaskiewicz et al. (1991), but it was rarely 
used in analysis of the observational data as well. 

Further, the cyclotron absorption that must take place 
near the light cylinder (Mikhailovsky et al. 1982) turns out 
to be so large that it will not allow the radio emission to 
escape the pulsar's magnetosphere (see, e.g., Fussel et al. 
2003). Finally, the limiting polarization effect had not been 
discussed seriously over many years, although it was quali- 
tatively clear that this effect must be decisive for explaining 
of high degree of circular polarization, typically 5-20%. 

Indeed, in the region of the radio emission generation lo- 
cated at 10-100 neutron star radii (these values result from 
the hollow cone model), the magnetic field is still strong 
enough so the polarization of the two orthogonal modes is 
indistinguishable from a linear one. For this reason, it is log- 
ical to conclude that the polarization characteristics cannot 
be formed precisely in the emission region, and the prop- 
agation effects are to play important role (cf. Mitra et al. 
2009). Nevertheless, in an overwhelming majority of the pa- 
pers, Eqn. (IT]) is used to investigate the polarization. 

Recall that the limiting polarization effect is related to 
the escape of radio emission from a region of dense plasma, 
where the propagation is well described in the geometrical 
optics approximation (in this case, the polarization ellipse 
is defined by the orientation of the external magnetic field 
in the picture plane), into the region of rarefied plasma, 
where the emission polarization becomes almost constant 
along the ray. This process was well studied (Zheleznyakov 
1977; Kravtsov & Orlov 1990) and was used successfully for 
numerous objects, for example, in connection with the prob- 
lems of solar radio emission (Zheleznyakov 1964). However, 
in the theory of pulsar radio emission, such problem has not 
been solved. Above the papers where the level r = resc at 
which the transition from the geometrical optics approxi- 
mation to the vacuum occurs, was only estimated (see, e.g., 
Cheng & Ruderman 1979; Barnard 1986), one can note only 
a few paper by Petrova & Lyubarskii (2000) (these authors 
considered the problem in the infinite magnetic field), by 
Petrova (2001, 2003, 2006), as well as the recent papers by 
Wang, Lai & Han (2010, 2011). 

The goal of our paper is to consider all three main effects 
(i.e., refraction, cyclotron absorption, and limiting polariza- 
tion) simultaneously in a consistent manner for realistic case. 
Not only the plasma density but also the magnetic field de- 
creases with increasing distance from the neutron star will 
be included into consideration. Also, the non-dipole mag- 
netic field, the drift motion of plasma particles, and realistic 
distribution function of outgoing plasma will be taken into 
account. 

In section 2 both ordinary and extraordinary waves 
propagation in the pulsar magnetosphere is briefly consid- 
ered. In addition, the hydrodynamic derivation of dielectric 
tensor of relativistic magnetized plasma is given. In section 
3 the main parameters of our model are discussed. In sec- 
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Figure 1. Dependence of the refractive indexes n on the angle 
between the wave vector k and external magnetic field B for 
Ap ^ 1. The lower branch corresponds to the O-mode. The angle 



tion 4 we discuss the transition from geometrical optics to 
vacuum propagation, the cyclotron absorption, and the wave 
refraction. The results of numerical calculations of the mean 
profiles of radio pulsar are presented in section 5. Finally, in 
section 6 we discuss our main results. 

It is necessary to stress that the main goal of this paper 
is in describing the theoretical ground of propagation effects 
only. For this reason, we are not going to discuss here in 
detail the theoretical predictions for real objects. This will 
be done in the separate paper. On the other hand, the theory 
described below is independent on the emission mechanism 
and, hence, can be applied to any theory of the pulsar radio 
emission. 



2 TWO ORTHOGONAL MODES 
2.1 On the number of outgoing Avaves 

As was already stressed, pulsar radio emission is highly 
polarized. The mean degree of linear polarization can reach 
40-60%, and even 100% in some subpulses (Lyne & Graham- 
Smith 1998). The analysis of the position angle demon- 
strates that in general the pulsar radio emission consists 
of two orthogonal modes, i.e., two modes in which their po- 
sition angles differ by 90°. It is logical to connect them with 
the ordinary (O-mode) and extraordinary (X-mode) waves 
propagating in magnetized plasma (Ginzburg 1961). 

Starting from the pioneering work by Barnard & Arons 
(1986), three waves propagating outward in the pulsar mag- 
netosphere were commonly considered (see, e.g., Usov 2006; 
Lyubarsky 2008). But in reality we have four waves prop- 
agating outwards. The point is that in the most of papers 
(see, e.g., Melrose & Gedalin 1999) the wave properties were 
considered in the comoving reference frame in which the 
plasma waves propagating outward and backward are iden- 
tical. But in the laboratory reference frame (in which the 
plasma moves with the velocity v ^ c) the latter wave is to 
propagate outward as well. 

As shown on Fig. ^ for zero angle between the wave 
vector k and external magnetic field B two of them, having 



refractive indices ni and 712, correspond to transverse waves. 
For infinite external magnetic field ni — n2 — 1. On the 
other hand, the waves ns and 714 corresponds to plasma 
waves propagating in different directions in the comoving 
reference frame. 

Moreover, as was demonstrated by BGI, it is the fourth 
wave 724 that is to be considered as the O-mode in the pulsar 
magnet osphere. It should be noted that it is valid for dense 
enough plasma in the radio generation domain for which 
Ap ^ 1, where 



< 7 > 



(2) 



Here and below cjp — {Ane^Ue/meY^'^ is the plasma fre- 
quency, Ue is the particle number density, rrie is the parti- 
cle mass, and 7 is the particle Lorentz factor. Further, in 
what follows we assume that the particle distribution func- 
tion Fg+ g_ (p) is one-dimensional. It results from the very 
high magnetic field in the vicinity of the neutron star where 
the synchrotron life time is negligible. In this case the brack- 
ets <> denote both the averaging over the one-dimensional 
particle distribution function and the summation over the 
types of particles: 



<(...)>= 



e I e 



(•••)^e+,e-(p)dp. 



(3) 



As a result, as shown in Fig.[T] for A^^ lit is the wave 
riA that propagates as transverse O-mode at large angles 
0^0^, i.e., at large distances from the neutron star. Here 



(jj'^j^ 



1/4 



(4) 



The second transverse wave for ^ ^ ^* is again the X-mode 
ni. Two other waves, 712 and ns, for which the refractive 
index n > 1, cannot escape from the magnetosphere as at 
large distances they propagate along the magnetic field lines 
(and due to Landau damping, see Barnard & Arons 1986). 
In the hydro dynamical limit one can easily obtain the 
dispersion curves shown in Fig. ^ from the well-known dis- 
persion equation in the limit of large magnetic field (see, 
e.g., Petrova & Lyubarskii 2000) 



fl — n^cos^6 



uj'^j^{l — nvcosO/c)'^ 



2 • 2 
-n sm ( 



0.(5) 



For <^ 0^ and for 0^0^ there are two transverse and two 
plasma waves, but for Ap ^ 1 the nontrivial transforma- 
tion from longitudinal to transverse wave takes place. This 
implies that in this case the mode 714 can be emitted as a 
plasma wave, but it will escape from the magnetosphere as 
a transverse one. 

As the refractive index 77,4 differs from unity, the ap- 
propriate ordinary mode deflects from the magnetic axis if 
^ ^ ^*. As was already mentioned, for the O-mode this ef- 
fect takes place if Ap > 1, i.e., for small enough distances 
from the neutron star r < ta, where 



^. ir>2 o A 1/3 1/3 T3I/3 

rA ^ 10 i?A/ 7100^12 ^ 



-2/3 



-1/3 



(6) 



Here R, P, and B12 are the neutron star radius, rotation pe- 
riod (in s), and magnetic field (in 10^^ G), respectively. Ac- 
cordingly, 7100 = 7/100, z^GHz is the wave frequency in GHz, 
and A4 = A/ 10^, where A = rie/riGj is the multiplicity of 
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A Observer 



Figure 2. X'Y' Z' frame connecting with the neutron star rota- 
tion (Z' axis) 



the particle creation near magnetic poles (noj = Q.B /2nce 
is the Goldreich- Julian number density). On the other hand, 
the transverse extraordinary wave with the refractive index 
ni ^ 1 (X-mode) is to propagate freely. As the radius ta is 
much smaller than the escape radius resc (Cheng & Ruder- 
man 1979; Andrianov & Beskin 2010) 



^. ir>3oA2/5 -6/5 c>2/5 -2/5 
resc -10 R\^ 7lOO ^12 ^GHz " 



-1/5 



(7) 



one can consider the effects of refraction and limiting polar- 
ization separately. In particular, this implies that one can 
consider the propagation of waves in the region r ^ resc as 
rectilinear. 



2.2 Extraordinary Avave 

Below for simplicity we assume that both two outgoing 
modes are generated at the same heights rem (few to tens 
NS radii), where the magnetic field can be considered as a 
rotating dipole 



B(<^) 



m(0) 3r 



■(m(<^),r) 



(8) 



Here = f^t is the corresponding pulsar rotation phase. 

In what follows it is convenient to use two coordinate 
systems. In a X'Y' Z' frame (Z^-axis is along fi; see Fig.^, 
we have 



m((/)) = sin acoscj) e^,/ + sin a sin e^/ + cosa e^ 



(9) 



In a XY Z frame (Z-axis is along the line of sight, Q, lies in 
XZ plane; see Fig. [sl we have 

m((/)) = (sin acosQcoscj) — sin ^cosa) e^c + sin a sin e^ 
+ (cosacosC + sinasinCcos(/))e2. (10) 

Therefore, in the XY Z reference frame the spherical angles 
'dm and (/)m of the vector m are (see Fig. Isl 




>Y 



Figure 3. XY Z frame connecting with the line of sight (Z axis). 
The dashed curve indicates the magnetic field line 



sm a sm c 



tan 0^ = — ^- : ^^— . (12 

sm Qcosa — sm q;cosCcos0 

In the rotating vector model (RVM) the p. a. is determined 
purely by the projection of magnetic field on the sky's plane, 
so it coincides with (j)m • The sign of the arctan term is deter- 
mined by the p. a. measured counter-clockwise in the picture 
plane, as is common in radio astronomy (Everett & Weisberg 
2001). 

As the aberration angle at the emission point is approx- 
imately ^rem/c, i.e., it is much smaller than the angular size 
of the emission cone I/7, we can easily find the position of 
the emission point, at which the magnetic field line is along 
the line of sight. This point rem = (^em, ^em, (/>em) in the 
XY Z frame is given by the spherical angles as 



^r, 



1 



- sini^ri 



^. 



(13) 
(14) 



Note that the impact angle (3 is the smallest angle between 
line of sight and magnetic moment m, is given by /3 = a — C- 
As a result, the trajectory of the extraordinary wave in the 
XY Z frame is given by the simple relation 



R = rem + TGz 



(15) 



cos'dyn = cosacosC + sin a sin ("cos(/ 



(11) 



This relation allows us to determine the magnetic field and 
all plasma characteristics along the ray. 

2.3 Ordinary Avave 

Let us briefly review the main points of the theory of 
ordinary wave propagation (Barnard & Arons 1986; Beskin 
et al. 1988; Petrova & Lyubarsky 1990b). In the geometrical 
optics limit, the equations of motion of a ray are 

dp± _ d k 
dl dk± rij ' 



(16) 
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dk± 
dl 



dpi 



(17) 



where p± is the distance from the magnetic dipole axis, / ~ r 
is the coordinate along the ray, the index ± corresponds 
to the components perpendicular to the dipole axis, e.g., 
0± — k±/k^ and Uj are the corresponding refraction indices. 
Below in this subsection for simplicity the plasma density 
is assumed to be independent on the transverse coordinate 

As the refraction of the 0-mode takes place at small 
distances from the neutron star r ^ r a (|6]) , the expressions 
for refraction indices can be borrowed from the theory in the 
infinite magnetic field when the dielectric tensor of plasma 
has a form: 










Here and below by definition 

cj = cj — (k, v). 



(18) 



(19) 



As the brackets <> denote the averaging over the particle 
distribution function, the singularities in the dielectric ten- 
sor coefficients in Cerenkov (and, below, in cyclotron) reso- 
nance vanish due to averaging over the wide particle energy 
distribution function. 

As a result, for the ordinary mode the equations take 
the following form: 



dy9± 
dl 


= 0^ 


^as 


2 


o± 




dO± 

dl 


_ 3 0± 
~ 4 


T 


as 



(as 



-f 



(5(^) + («B-^J^) 



1/2 



yo^B 



+ {aB-e^Y) 



1/2 



where as is the inclination angle of the magnetic field line 
to the magnetic axis. As was already mentioned, for large 
enough angles 0^0^^ the ordinary wave propagates 
rectilinearly as well. From this condition and the solution of 
the equation above one can find 



^±(00) = 



(^)* 



1 






/.0.36 
/em 



/r \0.15 

(if) -(22) 



Here cjpo is the plasma frequency near the star surface, and 
index 'em' corresponds to the quantities on the generation 
level. Besides, the dimensionless factor 



/ 



c 
OR 



R 



(23) 



where the angle Om is measured from the magnetic axis, 
determines the position of the radiation point within the 
polar cap. The angle 20 ^[00) then determines the angular 
width of the emission beam. Finally, the "tearing off' level 
It defined by the condition — 0^ is equals to 



2R 



It gives 
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IIR 

Figure 4. Model simulation of the refraction equation for the ab- 
sence (upper curve) and the presence (lower curve) of the trans- 
verse number density gradient 



Thus, this level locates much deeper than the level of the 
formation of the outgoing polarization resc ^ 1000 i? (I7|. 

As was already stressed, these results were obtained for 
the case of neglecting transverse plasma density gradients. 
To check the validity of this approximation, we show on 
Fig. |4] the characteristics of the wave propagation for two 
different density profiles. The upper curve corresponds to 
the absence of transverse gradients, and the lower one cor- 
responds to the case of the hollow cone distribution that 
will be used everywhere below. As they are quite similar, 
/ 20 the analytical results obtained above will be used. Also it is 
shown that the analytical estimates (22) and (24) are cor- 



rect enough. In more detail the procedure we have used is 
described in Appendix A. 



(21) 



2.4 Dielectric tensor 

For reasonable parameters of the plasma filling the pul- 
sar magnetosphere one can neglect the effect of curvature 
of magnetic field while considering the propagation of radio 
waves (see, e.g., Beskin 1999). On the other hand, as the 
level of the formation of the outgoing polarization resc vj 
locates in the vicinity of the light cylinder (Cheng & Rud- 
erman 1979; Barnard 1986), it is necessary to include into 
consideration the nonzero external electric field (Petrova & 
Lyubarskii 2000). In this paragraph z-axis is selected along 
the direction of magnetic field and the wave vector lies in 
xz-plane. 

Our goal is to find the permittivity tensor of relativis- 
tic plasma in perpendicular uniform magnetic and electrical 
fields. In the derivation we take into account the fact that in 
the strong enough magnetic field the unpertubated motion 
of particles is the sum of the motion along the magnetic field 
lines and the electrical drift in the perpendicular direction: 



M|b + U. 



(26) 



/. 



0.48 
em 



. APiTD dO.24 -0.48 -0.72 o0.24 a 0.24 p- 
40RP Z^GHz 7lOO ^12 A4 /e, 



/r \ -0.20 



/r \ -0.2 

(if) -(2^) 



Here b = B/B is the unit vector along the direction of 
magnetic field, and U = c[E, B]/B^ is the drift velocity. In 
what follows we will use another form of this equation 

Vo = [fi,r] + ci||B (27) 

resulting from the condition E + [V, B]/c = (BGI; Gruzi- 
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nov 2006). Here i\\ is the scalar function which will be de- 
termined below. 

To find the permittivity tensor Sij we have to find the 
motion of plasma particles in the homogeneous fields per- 
tubated by the plane wave. We start from linearized Euler 
equation: 



d_ 
dt 



+ Vo V Up = e Ue + 



^V 



,^B 






(28) 



(29) 



Writing now the clear relations for particle number density 
and electric current pertubations 

ddrie 



+ div(ne^V + ^rieVo) = 0, 



(31) 
(32) 



where aij is a conductivity tensor, one can determine the 
permittivity tensor by the following relationship (Ginzburg, 
1961) 



Sij + 



47rz 



(33) 



and the relation between fields in the electromagnetic wave: 



<5B 



-[k,5E]. 



(30) 



As a result, the expression for tensor Sij in the infinite mag- 
netic field looks like (the full expressions can be found in 
Appendix B): 



,2. .2,., 2^,2 



^P^U 






< 



< ^ > 

kzUyUJp{uJ-kxUx)lu 



a;2^3^,2 



>\ 



< 



kzUyUJp{uj-kxUx)lu 
-L~ <- ,~,2, ,.,3 > / 



(34) 



Here now 

uj = uj - kxUx - kzV\\ , (35) 

and 

7^ = (l-[/2/c')-'/'. (36) 

3 MAGNETOSPHERE MODEL 

3.1 Magnetic field structure 

As the formation of the outgoing polarization locates in 
the vicinity of the light cylinder, it is necessary to include 
into consideration the corrections to the dipole magnetic 
field which, actually, determines the disturbance of the S- 
shape form (fTj) of the p. a. swing. In this work we discuss the 
following models of magnetic field 



where here 



: Bd + B^ 



(37) 



where the field Bd connects with the dipole magnetic field 
of the neutron star, and the field Bw corresponds to the 
outgoing wind. 

For Bd we discuss two possible models. 

(i) The model of "non-rotating dipole", in which we ne- 
glect radiative corrections in the pre-exponential factors (see 
Landau & Lifshits 1975 for more detail): 



Br 



- sin sin a Re [E^ (r, (p,t)] + 



2^^cos^cosa, 



B. 



-^cos6'sinaRe[^x(r, (/?,t)] + 

|m| . , 

— — smt/cosa, 

— —^ sin a Re [iE^{r, (p,t)], 



(38) 



(39) 



(40) 



Qr 



E^{r, (p, t) = E^(r, ^p - Qt) = exp li \- i(p - iQt J . (41) 

As a result, since the ray propagates almost along the radius, 
i.e., r ~ ct, the full compensation of time and radial contri- 
butions in the factor exp (iQr/c -\- icp — iQt) takes place, so 
the ray does not feel the dipole rotation. 

(ii) The model of filled magnetosphere, which corresponds 
to a rigidly rotating dipole. 



Br 



- sin sin a Re [exp (icp — iQt)] + 



2^-^cos^cosa, 



-cos^ sin a Re [exp (icp — iQt)] + 



sin^cosa. 



■ sin a Re [i exp (icp — iQt)] . 



(42) 



(43) 



(44) 



Such a magnetic field at the distances r < i?L was obtained 
by BGI (1993) and by Mestel et al. (1999) as a consistent 
solution of the force-free equation describing neutron star 
magnetosphere. 

As to the wind component, we use the following ex- 
pressions corresponding to the so-called "split monopole" 
solution (Michel 1973; Bogovalov 1999) 

^tot 



B. 



B^ 



27rr2' 



^tot sin t 
r 



27ri?L 



(45) 
(46) 



Here ^tot = 27rfr{Q/c)\in\ is the total magnetic flux through 
the polar cap, and /r '^ 1 and f^^l are the dimensionless 
constants. We see that the former term describes the quasi- 
monopole radial magnetic field. Such a structure was ob- 
tained not only for the axisymmetric force- free (Contopoulos 
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Table 1. Models of the magnetic field structure 



model 


A 


B 


c 


dipole 


i 


ii 


ii 


wind 


U = o 


U = o 


U = i 



et al. 1999; Timokhin 2006) and MHD (Komissarov 2006) 
numerical simulations but it describes well enough the mag- 
netic field of the inclined rotator as well (Spitkovsky 2006). 
As we are actually interested in the disturbance of the dipole 
magnetic field inside the light cylinder only, we do not in- 
clude here into consideration the switching of the radial field 
in the current sheet in the equatorial region. As the total 
magnetic flux through the polar cap depends only weakly 
on the inclination angle a (BGI; Spitkovsky 2006), we put 
here for simplicity /r = 1 (for zero longitudinal current fr 
changes from 1.592 to 1.93). 

Besides, the latter term B^ (46) corresponds to the 



toroidal magnetic field connected with the longitudinal elec- 
tric current flowing in the magnet osphere. It is well-known 
that to support the MHD (in particular, force-free) outflow 
up to infinity the the total longitudinal current / is to be 
close to the Michel (1973) current /m = Q'^tot/4:7T (Con- 
topoulos et al. 1999). It corresponds to f^ ~ 1. On the 
other hand, for realizing this current for inclined rotator 
with dipole magnetic field, it is necessary to suppose that the 
current density j\\ is much larger than the local Goldreich- 
Julian current jgj ~ ^Bcosa/27v (Beskin 2010). As it is not 
clear whether the Michel current /m > Igj can be realized 
in the pulsar magnetosphere, in what follows the parameter 
f^ can be considered as a free one. 

In Table [l] we present the notation of the models which 
will be used in what follows. Magnetic field structure for 
model C (and for orthogonal rotator) is shown in Figis] It 
is qualitatively similar to the numerical model obtained by 
Spitkovsky (2006). 



3.2 Plasma number density 

Recall the well-known property of the one-photon par- 
ticle production in a strong magnetic field: the secondary 
particles are produced only if the photon moves at large 
enough angle to the magnetic field line (Sturrock 1971; Ru- 
derman & Sutherlend 1975; Arons & Scharlemann 1979). 
Since the relativistic particles near the neutron star surface 
can move only along the field lines with a Lorentz factor 
7 = (1 — v?i/c^)~'^^'^ {v\\ is the particle velocity along the 
magnetic field), the hard gamma-quanta emitted through 
curvature mechanism also begin to move along the field lines. 
As a result, the production of secondary particles will be su- 
pressed near the magnetic poles, where the magnetic field 
is nearly rectilinear. Therefore, one would expect the sec- 
ondary plasma density to be suppressed in the central region 
of the open field lines (see Fig. |6|. It is this property that 
lies in the ground of the hollow cone model. 

Bellow we assume that the plasma number density on a 
polar cap is known. It is convenient to rewrite it in the form 









(47) 
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Figure 5. Structure of the magnetic field lines in equatorial plane 
in corotating frame for the orthogonal rotator x = 90°- The mag- 
netic field is considered as the sum of dipole and "split-monopole" 
(model C). The scale corresponds to the light cylinder i^L = c/fl 




rJRo 



Figure 6. Plasma space distribution function g(r^) on the polar 
cap as a function of the distance r±_ from the magnetic axis for 
/o = 0.25. The dashed line corresponds to /o = 



Julian number density, i.e., it does not depend on the 
inclination angle a. Further, the multiplicity parameter 
A = ne/riQj determining the efficiency of the pair creation 
is (Daugherty & Harding 1982; Gurevich & Istomin 1985; 
Istomin & Sobyanin 2009; Medin & Lai 2010) 



A- lO'' 



10" 



(48) 



Here n^j = Q.B/2Tice is the amplitude of the Goldreich- 



Finally, the dimensionless factor giOm^^m) ^ 1 describes 
the real number density of the secondary plasma in the vicin- 
ity of the neutron star surface as a function of magnetic pole 
angles Om and Lprn- 

The procedure described below allows us to determine 
the properties of the outgoing radiation for arbitrary number 
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density Ue within the polar cap. For ihustration we consider 
an axiahy symmetric distribution 



9if) = 



P" + fl 



(49) 



where / = r\/ Rq is the dimensionless distance to the mag- 
netic axis, Rq — (Q.R/cY''^R is the polar cap radius, and 
the parameter /o describes the hole size in space plasma 
distribution (see Fig. p|. 

As was already stressed, in this paper we are going to 
discuss the theoretical ground of propagation effects only. 
For this reason we assume here for simplicity the intensity 
of radio emission in the emission region to be proportional 
to the number density of outgoing plasma. In reality the 
directivity pattern may differ drastically from the particle 
profile. 

To determine the number density rie of the outgoing 
plasma in the arbitrary point of the magnetosphere, we use 
quasi-stationary formalism, which is valid for quantities that 
are functions of 0— Qt. For such functions all time derivatives 
can be reduced to spatial derivatives by the following rules 
(Beskin 2009): 



l« =-"!«' 



id_ 

cdt 



V = Vx [/3r,V]-(VV)/3r, 



for any scalar (Q) and vector (V) functions. Here 

PR — . 

C 

Using now the continuity equation 
^''^+div(neVo) = 0, 



(50) 



(51) 



(52) 



(53) 



dt 
where the velocity Vo is given by Eqn. (27), one can obtain 

(BV)(nei||) = 0. (54) 

Hence, the product 

rieiw — const (55) 

remains constant along the field lines (BGI, Gruzinov 2006). 
Taking now into account only the first order by Q.r/c and 
assuming that the velocity of the outflowing particles is close 
to the light velocity c, we finally obtain 

1 



B 



[l-(b,/3R) 



(56) 



Thus, to determine the number density Ue at an arbi- 
trary point along the ray trajectory it is enough to know the 
number density and magnetic field B at the base of a given 
field line on the neutron star surface. But for this it is nec- 
essary to produce the back integration along the field line 
from any point along the trajectory up to the star surface. 

In Fig.[7|we show the dependences of the distance of the 
foot points to the magnetic axis r± within the polar cap as 
a function of the distance r along the ray for rotating dipole 
(model B) and for magnetic field including the monopole 
wind (model C) for the phase point (j) = —5°. As for model 
B the appropriate value can be obtained analytically (dashed 




Figure 7. Dependences of the distance of the foot points to the 
magnetic axis within the polar cap as a function of the distance 
r along the ray for rotating dipole (model B) and for magnetic 
field including the monopole wind (model C) for the phase point 
4> = —5°. Dashed curve corresponds to analytical expression (|57|> 



where iprn is the angle between local point on the ray vector 
and momentary magnetic axis, one can conclude that the 
precision of our procedure is high enough. 



3.3 Energy distribution and cyclotron absorption 

As was shown by Andrianov & Beskin (2010), for large 
enough shear of the magnetic field along the ray (as will be 
shown below, this condition does hold in the pulsar mag- 
netosphere), all the polarization characteristics of outgoing 
radiation depend on the diagonal components of dielectric 
tensor, that are not sensitive to the difference of the e^e~ 
distribution functions. This fundamental property allows us 
to consider the electron and positron energy distribution 
functions to be identical. Our particular choice is (see Fig. pi) 



f(7) 



670 7'' 

2i/67r 276+76- 



(58) 



This distribution has the maximum for 7 = 70 that is 
assumed as a typical Lor entz- factor of plasma particles, 
and has power-law spectrum 7"^ for large Lorentz-factors 
7 ^ 70. Thus, it models well enough the energy distribu- 
tion function obtained numerically (Daugherty & Harding 
1982; BGI). 

Finally, consider the cyclotron absorption taking place 
in the region where the condition cjb = 'yyuoj holds (see 
Appendix B for more detail). As is well-known, the cyclotron 
resonance locates at the distances 



9 1 n^ 7? 7.-1/3^-1/3 pi/3^-2/3 
2 1U itz^GHz 7ioo ^12 ^0.1 



(59) 



r± 



sinipri 



(57) 



comparable with the escape radius resc ^ (Mikhailovsky 
et al. 1982). It implies that these two effects are to be con- 
sidered simultaneously. On the other hand, as was already 
stressed, in this region one can neglect the wave refraction, 
i.e., to put Re [n] = 1. Remember that the estimate (59) was 



obtained for zero drift velocity [7 = 0. Nevertheless, as one 
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Figure 8. Particle energy distribution function 




Figure 9. Cyclotron resonance radius as a function of the pulsar 
rotation phase for u = IGHz, 70 = 100, and magnetic field model 
B 



can see on Fig. [9] the result of calculati on for real case is in 
qualitative agreement with estimation (59]F 



As a result, the intensity of outgoing radiation can be 
determined as 



/o exp( 



(60) 



where /o is the intensity in the emission region, and the 
optical depth r = 2(jj/c f Im [n] dl can be found using the 
clear relation 



Im [n] ^ Im [£yfyf]/2. 
As a result, we have 



(61) 



00 00 

j jn.{l)^F{^)5 



\^B\ 



-^ — '^uj I d7d/ 



ruec 



00 



\lobW\-U^Ic 



Al. 



(62) 



^ Here the delta-function for particle energy distribution was as- 
sumed, because in the case of distribution function ( |58| there is 
the wide zone of cyclotron resonance. 



Here we use the approximation vljc? ?^ 1 — U'^ /c? . In the 
case of identical distribution functions for electrons and 
positrons the difference in absorption of the O- and X- modes 
is proportional to rAN/N ^ t/A, that is almost negligible 
for A ::> 1. It should be noted that the same result was 
obtained by Wang et al (2010) using exact solution of the 
equations for the Stokes parameters near the cyclotron res- 
onance under the typical pulsar conditions. However, this 
result differs from one obtained by Petrova (2006). 

Remember that for evaluation one can use the simple 
relation (Mikhailovsky et al. 1982) 



T ^ A(l — cos^r. 



(63) 



Hence, for rres ~ 0.1 i?L, A ~ 10^, and ^res ~ 0.1 the optical 
depth is to be high enough (r ~ 10). On the other hand, as 
shown on Fig. |6] and Fig. |7| for rres ~ 0.1 i?L the ray passes 
the very central parts of the open field lines region where the 
plasma number density n© can be much smaller than Ancj 
idif) ^ !)• For this reason, as will be shown below, the 
absorption of the outgoing radiation can be not so strong. 

Finally, as it is shown on Fig. |9] in the case of rotating 
magnetosphere rres depends on the pulsar rotation phase. 
Competition of these two effects (i.e., non-uniform plasma 
density distribution and difference in cyclotron radii) deter- 
mines which part of the beam, leading or trailing one, will 
be absorbed more efficiently. 



4 LIMITING POLARIZATION 

The limiting polarization effect is well-known 
(Zheleznyakov 1977). When the radiation escapes into 
the region of rarefied plasma, the wave polarization ceases 
to depend on the orientation of the external magnetic 
field. At the same time, in the domain of the dense enough 
plasma where the geometrical optics approximation is 
valid, the orientation of the polarization ellipse is to be 
determined by the direction of the external magnetic field. 
This implies that the geometrical optics approximation 
under weak anisotropy conditions becomes inapplicable 
and the question about the pattern of the limiting polar- 
ization effect should be solved by using the equations that 
describe a linear interaction of waves in an inhomogeneous 
magnetoactive plasma. 

Traditionally to describe general radiative transfer in 
magnetoactive plasma four first-order differential equations 
(for all four Stokes parameters) are used (Sazonov 1969; 
Zheleznyakov 1996; Petrova & Lyubarskii 1990; Broderick 
& Blandford 2010; Wang et al. 2010; Shcherbakov & Huang 
2011). Budden eqution, i.e., the second-order equation to the 
complex function actually corresponds to the same approach 
(Budden 1972; Zheleznyakov 1977). On the other hand, both 
the standard and the Zheleznyakov-Budden approaches are 
not quite convenient for quantitative estimates of the polar- 
ization of the escaping emission in general case. But since 
we are going to describe the propagation of originally fully 
polarized waves, not the ensemble of waves, we actually need 
only two equations for observable parameters, i.e., the posi- 
tion angle and the Stokes parameter V. 

There exists a different approach that allows us imme- 
diately write down the equations for these observable quan- 
tities, namely, the Stokes parameter V, defining the circular 
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polarization and the position angle p. a., characterizing the 
orientation of polarization ellipse (Kravtsov & Orlov 1990). 
This approach is valid in the quasi-isotropic case, i.e., in the 
case when the dielectric tensor can be presented as 



£ij — edij + Xij 7 



(64) 



where the anisotropic part Xij is small as compared to 
isotropic one. In this case we have two small parameters 
— general WKB parameter 1/kL and 



An/ni,2 ^ Xij /ni,2 < 1- 



(65) 



As a result, the solution can be found by expansion over this 
two small parameters. 

As one can check, these conditions are just realized in 
the pulsar magnetosphere (Andrianov & Beskin 2010). In- 
deed, in the region r ^ resc "^ 10^ R the value of v = ^p/^^ 

v^lO-^X4Bi2iyGLP~\ (66) 

is much smaller than unity. Accordingly, the deviation of the 
refractive indices from unity, |ni,2 — 1| '^ 'y, is also very small 
here, so we can neglect the wave refraction in the polariza- 
tion formation region. 

The Kravtsov- Orlov equation 

^ + ^ [\Xba - Xab) + [Xba + Xab)COs2k^ - 



dl 



4c' 



iXaa -Xbb)sin2e], 



(67) 



is the equation for the complex angle G = Gi + iG2, where 
Gi is a position angle and G2 determines the circular polar- 
ization by the relation 



y = /tanh2G2. 



(68) 



Here / is the intensity of the wave. The components of the 
dielectric tensor Xij are to be written in a frame of unitary 
vectors a and b in the picture plane where a is determined 
by the projection of the vector Ve. Finally, 

K = l/2(a.[V,a] + b.[V,b]) (69) 

is the ray torsion (see Kravtsov & Orlov 1990 for more de- 
tail) . It can be easily understood that the rotation of position 
angle described by the ray torsion is fictions and describes 
only the rotation of coordinate system. As a result, we can 
write down 
dGi _ uj_. 

"dT ~ 2^' 

_lcj 
2V 
dG2 _ 1 cj 

"dT ~ 2 c^ 

Here / is a coordinate along the ray propagation, and the 
angle /3b (l) defines the orientation of the external magnetic 
field in the picture plane (defined as tan /3s = By/Bx, 
where Bx and By are the components of the magnetic field 
vector in a XYZ system, see Fig. |3|. Further, 



-Im[£3,/^/] 

^Acos[2Gi - 2/3b{1) - 2^(/)]sinh2G2, (70) 
-Asin[2Gi - 2/3s(/) - 2^(/)]cosh2G2. (71) 



A = =FW(Re[£^/^/])2 + 



(72) 



where the signs correspond to the regions before/after the 
cyclotron resonance and 

2Re[£a:'y'] 



tan{26) 



(73) 



Finally, s^/^/ are the components of plasma dielectric tensor 
in the frame where the z-axis directs along the wave prop- 
agation and the external magnetic field lies in the xz-plane 
(see Appendix C). As was already stressed, the singularities 



at the cyclotron resonance in equations (70)-(71) are absent 



due to averaging over wide particle energy distribution. For 



distribution function (58) this averaging can be done ana- 
lytically. 

We would like to note that in these equations the circu- 
lar polarization is defined as it is common in radio astronomy 
(positive V corresponds to LHC polarization) . Nonrelativis- 
tic version of the above equations is given in Czyz et al. 
(2007). It should be mentioned that the equation for the 
Stokes vector evolution has been recently shown to be de- 
rived directly from the Kravtsov-Orlov quasi-isotropic ap- 
proximation (Kravtsov & Bieg 2008). 

As one can see on Fig. |10| in the geometrical optics 
region Eqns. (70)-(71) describe oscillations of the angle Gi 



near the value Gi = /3b-\-S. As the ray moves into the region 
of rarefied plasma, the length of the spatial oscillations L ^ 
c/{(jjA n) increases and in the region r > resc becomes larger 
than the characteristic length r. As a result, the angles Gi 
and G2 become constant for r ^ resc- They are the values 
that characterize the outgoing radiation. 



Thus, the basic equations (70)-(71) generalize ones ob- 



tained by Andrianov & Beskin (2010) for zero drift velocity 
U = when Re[s^/^/] = and, hence, ^ = 0. In particu- 
lar, they now include into consideration the aberration ef- 
fect considered by Blaskiewicz et al. (1991). This effect was 
also considered by Petrova & Lyubarskii (2000), but for the 
infinite magnetic field only. It is important that in Eqns. 
(70)-(71) the angle Gi is measured relative to the labora- 



tory frame because these equations contain the difference 
between Gi and /Sb only. 

Equations above have the following important property. 
For homogeneous media {/3b = const, Sij = const) the pa- 
rameters of polarization ellipse Gi and G2 remain constant 
if the following conditions are valid: 



ei 



I3b+S, 



*-e. = '=^^-i ,74) 



Gi = Pb + S + 7r/2, sinh2G2 = 



Im [sx' 



y'\ 



A 



l_ 
Q 



Here (see the definition of e^/^/ in Appendix D) 



Q 



^^x' y' 



(75) 



(76) 



This closely corresponds to the polarization of the two nor- 
mal modes, the former corresponding to the 0-mode, and 
the latter to the X-mode. In addition, the following impor- 
tant property holds: irrespective of the pattern of change 
in plasma density and magnetic field along the trajectory, 
if two modes were orthogonally polarized in the beginning 
(G^^^ - of ^ = 7r/2, G^^^ = -G^^^), then this property will 
also be retain subsequently, including the region where the 
geometrical optics approximation breaks down. 

Finally, as was already mentioned, in the region r ^ resc 
one can put dGi/d/ ~ d(/3s +^)/d/. Hence, for high enough 
shear of the external magnetic field along the ray propaga- 
tion when the derivative d{pB + 5)/dx is high enough, the 
first term in the r.h.s. of Eqn. (70) may be neglected. As for 
G2 ^ 1 we have sinh2G2 ~ tanh2G2, one can write down 
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for V/I = tanh2e2 

y 1 d{PB + s )/dx 

l"^ \Q\ A 

Here 



1 



cos[2ei-2/3B(Z)-25(0]' 



I Ml 



l[/,/c)-cos6l(l-[/Vc')|. 



(77) 



(78) 



X = Q//c, and we used Xtigj for the plasma number density. 
Thus, the sign of the circular polarization will co- 
incide with the sign of the derivative d(/3B + 5)/dx for 
the 0-mode and they must be opposite for the X-mode. 
This approximation can be used for large enough deriva- 
tive d{l3 B -\- 5) / dx ^ 1 (i.e., for large enough total turn 
A (/3s + ^) ^ 1 within the light cylinder i?L = c/Q), and 
for small angle of propagation ^ ^ 1 through the relativis- 
tic plasma {v\\/c ^ 1). Both these conditions are valid in 
the magnetospheres of radio pulsars with a good accuracy. 
Indeed, assuming that U/c <^ 1 and Ux/c ^ U/c ^ one 
can obtain 



A', 



1 

2^ 



Ux 



U' 



2r 



<i. 



(79) 



So, the Stokes parameter V (77) is to be much larger than 



Vb = i^/Q resulting from standard evaluation (Ginzburg, 
1961). 

It is important that in this case in the region where the 
geometrical optics is valid the circular polarization is to be 



determined by the value of A ( 72 ) which does not depend on 



imaginary non-diagonal components of the dielectric tensor 
Xab and Xba- This fundamental property is well-known in 
plasma physics and crystal optics (see, e.g., Zheleznyakov et 
al. 1983; Czyz et al. 2007), but up to now it was not used in 
connection with the pulsar radio emission. For radio pulsars 
this property is especially important because for electron- 
positron plasma the imaginary part of the dielectric tensor 
depends significantly on the difference in particle energy dis- 
tributions which is not known with the enough accuracy. 

Finally, our numerical simulations show that the sign 
of the derivative d(/3s + 6)/dx is opposite to the sign of 
dp.a./d(j). As one can see from Eqn. (77), this results in an 
important prediction: 

• For the X-mode the signs of the circular polarization V 
and the derivative dp.a./d(/) should be the SAME. 

• For the 0-mode the signs of the circular polarization V 
and the derivative dp.a./dcj) should be OPPOSITE. 



This implies also that the effects of the particle drift motion, 
as was already found by Blaskiewicz et al. (1991) (see also 
Hibschman & Arons 2001), shifts the p. a. curve to the trail- 
ing part of the mean profile. As it is shown below, our results 
are in qualitative agreement with this statement. Moreover, 
at present there are some observational confirmations of this 
property (Mitra & Rankin, 2011). Nevertheless, the main 
distinction of our theory is in self-consistent definition of 
Tesc (and, i.e., p. a. shift value) on the direct solution of po- 
larization transfer equations, that depends not only on the 
geometry, but on plasma parameters as well. 

Typical evolution of angles Gi and G2 are presented on 
Fig. [To] It shows that the analytical estimate of the escape 
radius ([tJ is qualitatively correct. Nevertheless, it should 
be mentioned that it depends on the plasma multiplicity 
factor, that depends effectively on pulsar rotation phase due 
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Figure 10. Evolution of the angles 0i,2 for O-mode, 



-5°, 



A = 10^, 70 = 50, V = IGHz, and magnetic field model B. Their 
small oscillations in the region r <C Tesc can be seen for O2 only 
due to the different scales in the upper and lower panels. The dash 
lines correspond to the geometrical optics values Bi = fe +5 and 
62 = 1/2Q 



to non-uniform plasma number density distribution. It also 
shows that real values of G2 are indeed much larger than 
the corresponding standard value 1/2Q. 



5 RESULTS 

Thus, in this paper the arbitrary non-dipole magnetic 
field configuration, arbitrary number density profile within 
the polar cap, the drift motion of plasma particles, and their 
realistic energy distribution function are taken into account. 
It gives us the first opportunity to provide the quantitative 
comparison of the theoretical predictions with observational 
data. Using numerical integration we can now model the 
mean profiles of radio pulsars and, hence, evaluate the phys- 
ical parameters of the plasma flowing in the pulsar magne- 
tosphere. 

It is necessary to stress that the detailed discussion of 
the morphological properties of mean profiles resulting from 
different inclination and impact angles is beyond the scope of 
our consideration and are addressed to future papers. The 
goal of this paper is in quantitative analysis of the prop- 
agation effects on the polarization characteristics of radio 
pulsars. In particular, we try to determine how the plasma 
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Table 2. Parameters of the 'ordinary' pulsar 
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parameters affect the 5'-shape of the position angle swing 
and the properties of the mean profile. 



5.1 Ordinary pulsars 

At first, let us discuss the results obtained by numeri- 
cal integration of equations (70)-(71) for "ordinary" pulsar 
(its parameters are given in Table 217 Everywhere below the 
dashed curves on the intensity panel show the intensity pro- 
file without any absorption. As was already stressed, in this 
paper for simplicity we suppose that it repeats the particle 
number density profile shown in Fig. |6] If the dashed curve 
is not shown, then the absorption is fatal and only the origi- 
nal intensity (which is normalized to 100 in its maximum) is 
shown. The dashed curves on p. a. panels show the prediction 
of the RVM-model M. Finally, pulsar phase is measured 
in degrees everywhere below. 

On Fig. 
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we show the intensity loo (60) (left panel) 



and the p. a. swing (right panel) for extraordinary X-mode as 
a function of the pulsar phase (j) for "non-rotating dipole" 
without the wind component (model A); the drift effects 
are neglected as well. The circular polarization degree does 
not exceed one percent here and that is why this curve is 
not presented in this picture. It results from approximately 
constant Pb along the ray. For this reason, as we see, the 
p. a. curve is nicely fitting by the RVM model. 

Further, the upper solid line corresponds to /o = 0.25, 
and the lower one corresponds to /o = 0.0025. The lower 
intensity curve shows that for the very small core in the 
number density /o = 0.0025 {r±/Ro = 0.05) the absorption 
is fatal and the emission cannot escape from the magneto- 
sphere. As was already stressed, this property can be eas- 
ily explained. Indeed, for /o ^ 1 the rarefied region of the 
"hollow cone" is actually absent, and the rays pass the cy- 
clotron resonance in the region of rather dense plasma. On 
the other hand, for /o ~ 1 the number density in the region 
of the cyclotron resonance is low enough for rays to escape 
the magnetosphere without strong absorption. 

Below we consider model C as the main magnetic field 
model of our simulation. As to model B, one can find the 
examples of mean profiles basing on this model in Wang et 
al. (2010). The main differences between our results and ones 
presented in the paper mentioned above are in shifting of the 
p. a. curve to the trailing part (because of the particle drift 
motion) and in possibility of absorption of the leading part 
of the pulse (due to the non-zero toroidal magnetic field). 

On Fig. 
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we show the p. a. swing (lower panels), the 
intensity /oo (solid lines on the top panels), and the Stokes 
parameter V (dotted lines) as a function of the pulsar phase 
(j) for extraordinary wave for magnetic field model C and 
for various multiplicity parameter A. It is obvious that the 
absorption increases with increasing A. In most cases, the 
trailing part of the mean profile is absorbed (see Dyks et al. 
2010 as well). 

Besides, as one can see, the p. a. curves differ signifi- 



cantly from the RVM one (dashed lines) as A growing. This 
property can be easily understood as well. Indeed, as the 
escape radius resc ([7| increases as A^^^, for large enough A 
the polarization properties of the outgoing waves are to be 
formed in the vicinity of the light cylinder, i.e., in the region 
with quasi- homogeneous magnetic field (see Fig.jsj. 

Further, as was already mentioned, the drift effect 
causes the p. a. curve to be shifted to the trailing part of the 
mean profile. It is necessary to stress that the opposite shift 
is to take place if we neglect the drift effect on the dielectric 
tensor (Andrianov & Beskin 2010; Wang et al. 2010). One 
can note that for high values of multiplicity parameter (i.e., 
for full absorption of the trailing part of the mean pulse) the 
observer will detect approximately constant p. a. Finally, as 



one can see from Eqn. (78), the maximum of circular polar- 



ization is also shifted to the trailing side, as larger devia- 
tions from the 5'-shape produce larger circular polarization. 
It is not visible on this picture because the trailing side of 
the beam is absorbed. Thus, one can conclude that the self- 
consistent quantitative analysis of observational data is to 
include these effects into consideration. 

Another point should be mentioned. It is known that 
correct accounting of the inverse Compton scattering of the 
surface X-rays can reduce the pair multiplicity significantly 
(see, e.g., Hibschman & Arons 2001). In this case plasmas 
effects on emission polarization can be negligible. 

On Fig. [13] we show the same dependences for various 
Lorentz- factors of outgoing plasma 70 = 10, 50, and 300. As 
the escape radius resc (M) decreases as 7~^^^, the largest shift 
of the p. a. curve takes place for small 70 = 10. Finally, on 
Fig. [M] one can see the same dependences for various wave 
frequencies u = 0.03, 0.2, and 0.5 GHz. As resc oc v~^l^ ^ the 
largest shift of the p. a. curve takes place for small frequen- 
cies. One can note that the full investigation of frequency 
dependence of mean profiles of radio pulsars is more com- 
plicate and must include, e.g., the detailed analysis of fre- 
quency dependence of the emission radius rem (BGI). This 
is beyond the scope of the article. 

As was demonstrated above, in general the sign of the 
circular polarization remains constant for a given mode. But 
under certain conditions the change of the V sign in the same 
mode may occur. It can take place when we cross the direc- 
tivity pattern in the very vicinity of the magnetic axis. Such 
an example for the X-mode is shown on Fig. [15] for P = Is, 
( = 49°, a = 48.5°, A = 10^ 70 = 50, and rem = 100 i?. 
Certainly, this example is illustrative only. To explain the 
change of the sign of the circular polarization in the centre 
of the mean profile the separate consideration is necessary. 

Finally, in Table l3] we present the values of maximum 
derivative (dp.a./d(/))max, that is commonly used for deter- 
mination of the inclination angle a (see, e.g., Kuzmin & 
Dagkesamanskaya 1983; Malov 1990; Everett & Weisberg 
2001). As we see, these values significantly depend on the 
plasma parameters (and can differ drastically from the RVM 
value). Hence, more precise specifying of magnetospheric 
and plasma model is necessary for quantitative analysis of 
the pulsar characteristics. 

In our opinion, this technique can be applied in future 
for radio pulsars for which the whole S'- swing is detected. 
Otherwise, the situation like in Fig. 



12 



for A = 10 may 
occur, when one can detect only constant part of the p.a. 
swing. At any way, one can conclude that if one try to fit 
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Figure 11. The intensity loo ( |60| (left panel) and the p. a. swing (right panel) as functions of the pulsar phase cf) (in degrees) for "non- 
rotating dipole" (model A). Here and below the dashed curves on the intensity panel show the intensity profile without any absorption. 
The upper solid line corresponds to /o = 0.25, and the lower one - to fo = 0.0025 



the whole p. a. curve by the RVM function (IT]), then for low 
enough number density (A < 10^) and high enough particle 
energy (70 > 50) the angles a and /3 obtained will satisfy the 
reality. But in other cases unexpected result may occur. As 
shown on Fig. 16 left panel (see Table [2] as well), the sim- 
ulated curve could be nicely fitted by the RVM model with 
non-zero shift value (fl]), but with the angles /3 ~ — 5° and 
a ^ 97° that drastically differ from real ones a = 45° and 
/3 — —3°. On the other hand, the precision in determination 
of the angles a and /S may be not so high (see Fig. 
panel) . 
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Table 4. Parameters of the millisecond pulsar 



p 


Bo 


a 


/3 


/o 


''"em 


70 


A 


20 ms 


108 G 


45° 


-3° 


0.04 


1.25R 


50 


10^ 



field. As a result, the cyclotron resonance takes place in the 
region of rarefired plasma for the trailing part of the pulse. 
Also, stronger deviations from the 5'-shape oi p.a. swing are 
found as compared to the case of ordinary pulsar, because 
the polarization forms closer to the light cylinder. 



5.2 Two modes profiles 

On Fig. [T7| an examples of the mean profiles including 
two orthogonal modes are presented. The ratio of the inten- 
sities in the radiation domain r = rem was assumed to be 



1/3. The jumps in the p. a. curves were done at 



r(0) /r(0) 

the phase where Iq = /x- As we see, this jump can dif- 
fer from 90°. This results from the different trajectories of 
two orthogonal modes. It is necessary to stress that in our 
consideration the ratio /q //x i^ ^^^^' ^^^ i^^ precise de- 
termination is a question for the radio emission generation 
theory. 

We see that for radio pulsars for which two modes can 
be observed the mean profiles are indeed to have the triple 
form, which, in general, are not to be symmetric. As the 
0-mode deviates from magnetic axis, we have to see the O- 
mode in the leading part, the X-mode in the centre, and 
again the 0-mode in the trailing part of the mean profile. 
The detailed comparison with observational data will be pre- 
pared in the separate paper. 



5.3 Millisecond pulsars 

As the last example, on Fig. [18] we show the mean 
profiles obtained for millisecond pulsar (the parameters are 
given in Table [4|. One important feature appearing here is 
that the leading, not trailing part of the mean profile can be 
absorbed. It is caused by the bending of the open field lines 
tube near the light cylinder due to non-dipolar magnetic 



6 DISCUSSION AND CONCLUSIONS 

In this paper we study the influence of the propagation 
effects on the mean profiles of radio pulsars. The Kravtsov- 
Orlov approach allows us firstly to include into considera- 
tion the transition from geometrical optics to vacuum prop- 
agation, the cyclotron absorption, and the wave refraction 
simultaneously. Arbitrary non-dipole magnetic field config- 
uration, drift motion of plasma particles, and their realistic 
energy distribution were taken into account. Using numer- 
ical integration, we found how the propagation effects can 
correct the main characteristics of the hollow cone model. 

To summarize, one can formulate the following results: 

• We confirm the one-to-one correlation between the signs 
of circular polarization V and the position angle derivative 
(dp.a./d(/)) along the profile for both ordinary and extraordi- 
nary waves. For the X-mode the signs should be the SAME, 
and OPPOSITE for the 0-mode. 

• The standard 5'-shape form of the p. a. swing (1) can 
be realized for small enough multiplicity A and large enough 
bulk Lorentz-factor 7 only. In other cases the significant dif- 
ferences can take place. The location of the p. a. maximum 
derivative (dp.a./d(/))max is shifted to the trailing side of the 
pulse. 

• The value of p. a. maximum derivative (dp.a./d(/))max, 
that is often used for determination the angle between mag- 
netic dipole and rotation axis, depends on the plasma pa- 
rameters [and differs from rotation vector model (RVM) 
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Figure 12. The same for model C and for different multiplicity factors A = 10^, 10^, 10^, and 10^. Here 70 = 50 and v = IGHz. Dotted 
lines correspond to Stokes parameter V 
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Figure 13. The same for various Lorentz-factors 70 = 10, 50, and 300. Here A = 10*^, and v = IGHz 



value] and, hence, cannot be used without more precise spec- 
ifying of magnetospheric plasma model. 

• In general, the trailing side of the emission beam is 
absorbed. 

In our opinion, even the preliminary consideration 
demonstrates that these conclusions are not in contradic- 
tion with observational data. The first statement is in good 
agreement with observations (Han et al. 1998; Andrianov & 
Beskin 2010). In some cases, the sign reversal in the core can 



occur (see Fig. 15 and Wang et al. 2010), that is detected 
for several pulsars (see also Radhakrishnan & Rankin 1990). 
The second point is in agreement with empirical model of 
Blaskiewicz et al. (1991). On the other hand, as the shift 
value depends on plasma parameters (assuming fixed geom- 
etry), the solution of the inverse problem (i.e., the fitting of 



the real pulsar profiles) can in principle provide us a chance 
of correct enough estimating of plasma parameters. Finally, 
the last point is in agreement with observations as well (see, 
e.g.. Backus et el., 2010). In more detail, it takes place for 
large enough multiplicity parameter A > 10^. But in some 
cases the leading part can be absorbed as well. It happens 
when the polarization forms close to the light cylinder. 

Finally, Fig. ^] demonstrates our understanding of the 
X- and 0-modes propagation in the pulsar magnet osphere. 
Here the dot-dash lines going to the panels with pulsar pro- 
files show different intersections of directivity pattern, that 
form the observable profiles. Resulting from the different 
propagation, the O- and X-modes produce two concentric 
cones, the inner (core) part corresponding to the X-mode 
and the outer (conal) one to the 0-mode. Remember that 
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Figure 14. The same for model C and for various frequencies v = 0.03, 0.2, and 0.5 GHz. Here A = 10^, and 70 = 50 



the analysis of the relative intensity in O- and X-modes is a 
question for generation theory. In the context of the present 
propagation theory this parameter should be considered as 
a free one. 

Another source of uncertainty which can affect the for- 
mation of the mean profile (and, in particular, results in 
the formation of its more complicated shape) is the possi- 
ble "spotty" of the directivity pattern, i.e., the presence of 
the separate radiative domains within the diagram (Rankin 
1993). As is well-known, such separate domains are observed 
as a subpulse drift (Ruderman & Sutherland 1975, Desh- 
pande & Rankin 2001). 

In addition, the full analysis is to include the possible 
linear depolarization. E.g., depolarization at the edges of 
the profile can be easily interpreted as a mixing of X and 



0-modes (Rankin & Ramachandran 2003). But the total 
degree of linear polarization (e.g., if we clearly see only one 
mode) is a question to propagation in the generation region. 
If this region is rather thick then one can expect the possible 
depolarization. For this reason we do not treat here the total 
degree of linear polarization at all. 

As we see, in general if only one mode dominates in 
pulsar profile, then we expect statistically D (double) mean 
profile for the 0-mode and S (single) one for the X-mode. Of 
course, we suppose here that there is nonzero natural width 
of the radiative domains that efficiently smooths the hollow 
cone corresponding to the narrow X-mode. Nevertheless, as 
is shown on Fig. [l9] the X-mode can in principal produce D 
profile and 0-mode can produce S ones. The Xd and Xs 
profiles can be obtained when only X-mode is produced in 
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Figure 15. The results of simulations for X-mode, P = Is, /S = -0.5°, a = 48.5°, A = 10^, 70 = 50, and rem = 1007?. It shows that 
under certain conditions the change of V sign in the same mode may occur 
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Figure 16. The left panel shows the p. a. swing (points) for 70 = 50, A = 10^, v = IGHz, a = 45°, and /3 = —3°. Solid line corresponds 
to the best RVM fit giving unrealistic values a = 97° and (3 = —5°. The right panel shows x^ map for RVM fitting 

Table 3. The position angle maximum derivative (dp. a. /d(/)) max and the shift Acf) of its position 



A 


70 


u (GHz) 


(dp.a./d(/))max 


RVM 


Act>n 


102 


50 




-9.47 


-10.14 


4.7 


104 


50 




-14.47 


-10.14 


7.4 


10^ 


50 




-11.70 


-10.14 


10.5 


103 


10 




-12.72 


-10.14 


11.5 


103 


50 




-9.86 


-10.14 


4.3 


103 


100 




-9.47 


-10.14 


4.7 


103 


300 




-9.46 


-10.14 


4.0 


103 


50 


0.03 


-15.92 


-10.14 


8.4 


103 


50 


0.5 


-11.82 


-10.14 


3.8 


103 


50 


0.2 


-18.02 


-10.14 


5.5 



18 V. S. Beskin and A. A. Philippov 



/''~'--,^-"'"\ 


/ \A=10^ 7-50 


/ \ 'o 


/ \ 


/ \ 


f \ 


1 \ 


/ \ 


1 \ 


/ \ 


/ \ 


y — --_'' ^_-- — -s^ 


/ \ 


/ \ 






// \ \ 


// V \ 


V ^^--^ \ 


// \. \ 


// \ N 


// \ \ 


/ / X^ \ 







Figure 17. Main profiles including both modes (O and X ones) for A = 10*^ and 10^. Here 70 = 50, and u = IGHz 



the formation region. On the other hand, if there are two 
modes with the similar intensity, the triple (T) and mul- 
tiple (M) profiles can be obtained as well. So, this picture 
is in agreement with Rankin's morphological classification 
(Rankin 1983). Our investigation provides an additional in- 
formation to morphological structure, i.e., the information 
of the type of the mode that forms mainly the mean pro- 
file. As was demonstrated, only the correlation between p. a. 
derivative and V signs provides us this chance. 

Finally, the determination of the angular dimension of 
the core and conal beams depends greatly on the emission 
radius (i.e., again, on the generation mechanism). So the de- 
tailed analysis is also beyond the scope of current investiga- 
tion (one can find empirical results on this in, e.g., Maciesiak 
& Gil 2011). 

Thus, the approach developed in this paper allows the 
predictions of the theory of radio emission to be quanti- 
tatively compared with the observational data. Moreover, 
based on the observational data, one can reject a number 
of models in which the shapes of the profiles of the position 
angle and the degree of circular polarization never realized 
in practice are obtained. Thus, it becomes possible to solve 
the inverse problem, i.e., to determine the parameters of the 
outflowing plasma and the magnetic field structure. 

We will be glad to any collaboration with observers for 
investigating the observational features of the mean profiles 
on the base of our theory. The source code of our calculations 
(in Mathematica 7 or C) is available under request. 
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APPENDIX A: REFRACTION 

In this work we consider the simple model of the refrac- 
tion obtained under the assumption that the plasma number 
density is constant within the polar cap, i.e., g{Om, ^m) — 1- 
To include refraction into consideration we introduce the 
"imaginary source" of radiation giving the same trajectory 
at large distances whence emitted parallel to the magnetic 
field line (see Fig. Al). As the "tearing off" level locates 
deeply in the magnet osphere, i.e., ta ^ i?L, one can use the 
analytical expression (22) for the angle ^±00 (BGI 1993). It 



gives for the polar angle of the emission point 



R 



0±c 



2 



(Al) 



(for dipole magnetic field it does not depend on the radial 
distance). As a result, we obtain for the polar angle of imag- 
inary source 
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Figure 18. The results of simulations for X-mode for the millisecond pulsar for multiplicity parameters A = 10"^ (left panel) and A = 10^ 
(rigth panel) 




Figure Al. Imaginary source corresponding to rectilinear prop- 
agation of the wave 



After some algebraic calculations, one can find the trajectory 

R = m((/))pcos^im + m ±{(l))p sin 0[m + rcz. 
Here 

\sm/3 — cos/3 tan ^^o 



pW 



I tan^im — tan^^< 
is the radius of imaginary source, and 



^ ^ 1 . /sin^ 

/3=- + -arcsmf^ 



Vl + tan^ 6'in 



where 
^ = 0i 






1/4 



(A3) 



(A4) 



(A5) 



Finally, in^((/)) is a unit vector perpendicular to in((/)) lying 
for every pulsar phase (j) in the plane containing the magnetic 
m((/)) and the wave k vectors. 



tan^ 



tan^^c 



3l + y/l + 8/9tan2 6>±^ 



APPENDIX B: DERIVATION OF DIELECTRIC 
TENSOR 



Since all the quantities in equation ( 28 ) are proportional 



(A2) to exp {—iujt + ikr), its solution can be easily obtained: 



6vx 



^B - 7^7^^^ 



-iuj'y ( 1 + 7^ ^ j Di + Iujb+ iCo^lu "" ^ ^ ) ^2 



(Bl) 
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Figure 19. Formation of the directivity pattern in the pulsar magnet osphere. The X-mode propagating rectihnearly forms the core 
component, and deflected from the magnetic axis O-mode forms the conal one. The lower curves on each top panels indicate the Stokes 
parameter V, and the bottom panels shows the p. a. swing. Corresponding pulsar proflles was taken from Weltevrede & Johnston (2008) 
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D2 = 
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Prom (|3T]) one can find the number density pertubation After making substitution to l|32| we obtain for dielectric 



Sue = -r- (k, (5v) . 



(B6) 



tensor components: 



^2^3^2 ^2(^^^ _ ^2^2^^2^ 



(B7) 
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Here the quantity cjq = cj — kzV\\ . Other components of 
the permittivity tensor are defined by its hermiticity. There- 
fore, in the expression for gyrofrequency 



iOB 



eB_ 

rrieC 



we should take into the account the sign of the charge e. One 



can check that moving to the appropriate reference frame 
this dielectric tesor reduces to the well-known one (Godfray 
et al. 1975; Suvorov & Chugunov 1975; Hardee & Rose 1975) 
obtained for the case E = 0. Finally, as we see, the cyclotron 
resonance corresponds to the condition ujb — 77t7CJ. 

As a result, one can find for the coefficients in Kravtsov- 



Orlov equations ( 70 )-( 71 ) 



A = 

^^ \px'y'\ — 
' y' ^x'x' 

tan(2^) = 

tan(^) = 



"2 ^ /r,2^/3 ,,,2 



2 2 



UJ^^'' (jO^ — j'^^IjUJ'^ 



> [(sin 0-U:,/cf + cos^O Uy/c^ 



< - 



2 2 



CJ^^S ^'^ _ ^2^2^^2 



> Uy/c (sin 6* - Ux/c) cosO, 



< - 



2 2 



LO'^J^ CJ^ 



■ 7^7^^^ 



>[(s 



■ Ux/c) — COS OUy/c 



lU^B 



cj cj| - 7^7^^^ 



> [cos^7^ — v\\/c{l — s\nOUx/c)\ 



2Uy/ccosO{smO - Ux/c) 
' (sin 6* - Ux/c)^ - cos^ OU^/c^ ' 

cos^ Uy/c 
sinO — Ux/c' 



(B14) 
(B15) 
(B16) 
(B17) 
(B18) 
(B19) 



APPENDIX C: ROTATION OF COORDINATE 

SYSTEMS 

The equations of the transition of the dielectric tensor 
from the coordinate system formulated in Sect. 3 to that 
used in Sect. 4 are the following 



■^y'y 



£y> 

^x' y 
£y'z 

^x' Z' 



Scccccos + Szz sin — {sxz + Szx) sin cosO, 



Sxx sin + Szzcos + (exz + £zx) sin cosO, 



uCOsO 



, sinO, 



= Syx sin + SyzCosO, 



= (Sx 



) sin cosO + SxzCos — Sz 



: sin 9. 



Here is the angle between the ray propagation and local 
magnetic field directions. Other components of the tensor 
can be defined by its hermiticity. 



APPENDIX D: NATURAL MODES 

For obtaining naturals modes of tensor e^/^/ it is con- 
venient to rotate the coordinate frame in xy-p\a,ne to make 
the condition Re [Sx'y'] = to be valid. One can find that it 
can be done by rotation defined by the following matrix: 



M = 



cosa 


— sin a 





sin a 


cosa 











1 



where 



tan(2a) = tan(2^) 



2Re [€x 



(Dl) 



(D2) 



In this frame the corresponding expressions for tensor com- 
ponents (defining as Ci/j/) are: 

Sx'x'cos a + £^/^/ sin a + 2Re [s^^/^/] sinacosa, 
y'y' — £3^/3;/ sin a -\- £yfyfCos fl " 2 Rc [s 3^ / y / ] slu fl cosfl , 

x'y 



e 

Cyf 

^x' 



ilYn[exy]- 



The naturals modes for this tensor (with pure imaginary 
x^y^ component) are well-known (Andrianov & Beskin 2010; 
Wanget al. 2010): 



Gi = /3, sinh2e2 = - 

Gi = /3 + 7r/2, sinh2G2 = 
where by the definition 



1 

Q' 
1 



^^x' y' 

It can be easily shown that in our case 
1 _ Im[e^/j,/] 



Q 



A 



(D3) 
(D4) 

(D5) 
(D6) 



Hence, Eqns. ( 74 )-( 75 ) really define the natural modes in the 



laboratory frame. The simple expression, which corresponds 
to a zero plasma temperature for Q is the following (cf. 
Melrose & Luo 2004 for zero drift): 



Q 



XujBUJ [(sin 6* - Ux/c)^ + U^/c^cos^e] 
273cj2 [cos6'(l - U^/c^) - vi\/c{l - smOUx/c)] 



(D7) 
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